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ABSTRACT 

Gradients of the sea surface topography (SST) -i.e., heights of the sea surface in relation 
to the geoid - are essential for the real-time modelling of ocean dynamics. Ocean current 
measurements indicate the existence of SST gradients as large as ±0.1 m per 10 2 km. A pre- 
requisite for remote sensing SST from satellites is a geoid model with at least ±6 cm resolu- 
tion through equivalent wavelengths. 

The only potential source of such data is a satellite-determined gravity field model. The 
internal statistics of the best such model available at present (GEM 9) indicate that favour- 
able signal-to-noise exists for the recovery of the dominant parameters of the quasi-stationary 
dynamic sea surface topography from GEOS-3 altimetry. 

The 1977 altimetry data bank available at Goddard is analyzed for the geometrical shape 
of the sea surface expressed as surface spherical harmonics after referral to the higher refer- 
ence model defined by GEM 9. The resulting determination is expressed as quasi-stationary 
dynamic SST. Solutions are obtained from different sets of long arcs in the GEOS-3 altim- 
eter data bank as well as from sub-sets related to the September 1975 and March 1976 equi- 
noxes assembled with a view to minimizing seasonal effects. 

The results obtained are compared with equivalent parameters obtained from the hydro- 
static analysis of sporadic temperature, pressure and salinity measurements of the oceans and 
the known major steady state current systems with comparable wavelengths. 
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The most clearly defined parameter (the zonal harmonic of degree 2) is obtained with 
an uncertainty of ±6cm. The preferred numerical value obtained (-43 cm) is smaller than the 
oceanographic value (-46) largely due to the effect of the correction for the permanent Earth 
tide. Similar precision is achieved for the zonal harmonic of degree 3. The precision obtained 
for the fourth degree zonal harmonic reflects more closely the accuracy expected from the 
level of noise in the orbital solutions, being a factor of 3 inferior to the values quoted above. 

Attempts to obtain the harmonics f slll and f sll0 were not successful because of the 
masking effect of ..he non-geocentricity of the system of reference used. The dominant ef- 
fect is a southward displacement of 1 .5 m along the polar axis. 

The results presented in this paper are preliminary. While some further progress of a 
limited nature may be forthcoming with improvements in the definition of orbits, the most 
important requirement for significant advances in remote sensing surface ocean dynamics 
using altimeter data, is the refinement of low degree tesseral harmonics of the satellite- 
determined gravity field model to 2 parts in 10 9 . 
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DETERMINATION OF SOME DOMINANT PARAMETERS 
OF THE GLOBAL DYNAMIC SEA SURFACE 
TOPOGRAPHY FROM GEOS-3 ALTIMETRY 


1. INTRODUCTION 


The basic equations of non-tidal motion of the surface layer of the oceans, excluding 
wave motion, can be represented on anx, x 2 coordinate system in the local horizon (x t 
oriented east, x 2 oriented north) by (e.g., see Mather 1976 for a review of the derivation) 


..re. d£« 1 3p, 

x ‘ + « k i jXk = ■ g ^‘^a^ +Fi 

obtained when i takes the values 1 and 2, f being the Coriolis parameter given by 

f = 2 co sin <t> 


6 y is the Kronecker delta, 


e ijk 


1 if 
-1 if 
0 if 


subscripts are ordered 1231231. . . 
subscripts are ordered 132132. . . . 
i=j; i=k, j=k 


( 1 ), 


( 2 ), 


(3), 


x*j, Xj, Fj being components of accelerations, velocities and frictional forces acting on the 
surface layer of the oceans in the directions x i? g is local gravity, p a the atmospheric pressure 
on the sea surface, p w the density of sea water and the dynamic sea surface topography 
(SST). 


f s is defined as the radial departure of the sea surface from the geoid. As the term “sea 
surface topography” has been used in the past to refer to the height of the sea surface above 
any arbitrary reference surface, it has been considered necessary to qualify by the use of 
the additional adjective “dynamic” as in the title of this paper. However, the abbreviation 
“SST” will always be used in this paper to refer to the quantity f s , as defined above. The 
height of the sea surface above the adopted reference surfaces will be called the height 
anomaly, designated by f if S is a rotating equipotential ellipsoid of revolution, or f ' if S 
were the higher reference model (Mather 1974, p. 90 et seq.). 

The currents in the surface layer of the oceans can be as large as 2 x 10- cm s -1 . Table 
1 summarizes the factors necessary to maintain a current of 10 cm s -1 at latitude 45°. Figure 
1 is a representation of the dynamic sea surface topography as produced from steric anom- 
alies in the oceans, compiled primarily from (Wyrtki 1975. Reid et al. 1977, and Usitzin 1974) 
Such global representations are not instantaneous, being based on sporadic measurements of 
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Table 1 


Magnitudes (Rounded Off) of Factors Which Can Maintain a Surface Layer 
Current Velocity of 10cm s" 1 at Latitude 45° 


Factor 

Required Magnitude 

Sea Surface Topography Gradient ,/dx,, 

0.4arcsec (20cm per 10 s km) 

Atmospheric Pressure Gradient dpj/dx*, 

30 mb per tO 2 km 

Wind Stress* 

40 ms' 1 (83mph) 


*F a computed Uom the relation 


H * De.-rth of mixed laye' (K^m) 


temperature and salinity in the oceans. It is estimated that the contours can be in error by 
up to ±20 cm due to temporal variations. This model is merely a useful guide in the design 
of techniques to recover the SST from satellite altimetry. It should also be noted that the 
contours shown contain a zero degree term of + I.l4m due to the oceanographic datum 
which is of no geodetic consequence. 

Fast fiowing quasi-st at ionary mid-latitude currents like the Gulf Stream in the western 
No r th Atlantic (e.g., NOAA 1975) and the Kuroshio in the western North Pacific (e.g., 
Cheney 1977) are maintained by steep sea surface topography gradients in excess of 1 m per 
10 2 km but with relatively short wavelengths. Other more moderate current systems are 
maintained by the interplay of wind, sea surface topography and frictional forces acting on 
the surface layer due to continental shelf margin bottom topography. The effect of SST on 
ocean circulation is a function of latitude (equation 1) being a minimum in equatorial waters, 

The GFOS-3 altimeter had a design criterion of ±50cm resolution in the short pulse 
mode. The comparison of overlapping passes of such altimetry (Mather and Coleman 1977, 
Table 2) showed that the altimeter appeared to provide a resolution of at least ±20 cm on a 
relative basis, for features with wavelengths greater than 30km, assuming that the discrep- 
ancies were not of oceanographic origin. Spectral analysis of pairs of overlapping passes in- 
dicated that non-trivial strengths of signal were obtained for several wavelengths in excess of 
100km (Mather 1977, p. 25). Further details of investigations of this type using a wider 
data base, are reported in (Mather et al. 1978b). 

These studies indicated that a basis existed for the recovery of limited information 
about the global SST from the short pulse mode G EOS-3 altimetry, despite the following 
shortcomings in the data: 












MERCATOR PROJECTION 


finite I Quasr-ytauonary Dynamic Sea Surface Topography Rela..*e to lOOOdb Surface 
g Wavelengths > 10 3 km 




(i) The orbits (and hence the values of t provided by Wallops Flight Center) were 
subject to error. These were as much as two orders of magnitude worse than the 
relative precision of the altimetry data, especially for the data collected in 1975. 

In 1976, the orbital errors were about one order of magnitude worse than the 
resolution of the altimeter. 

(U) The GEOS-3 altimetry data was collected as a set of finite passes with lengths up 
to 9000 km. No attempt was made to obtain near-simultaneous global coverage 
as envisaged in the case of SEASAT-A. 

A study of overlapping passes in the Tasman and Coral Seas (Mather et al. 1977, p. 36) 
showed that the root mean square (rms) discrepancy between sea surface heights on two 
passes of length greater than 3000km dropped on the average to around ±30cm or less when 
the radial orbital error was modelled by a bias (b) and a tilt (c). This practice has been widely 
used by more than one GEOS-3 Kincipal Investigator when attempting to refine the altim- 
eter data for various srudies. 

The monthly analyses of GEOS-3 altimetry data in the Sargasso Sea (Mather et al. 
1978b) shows i. models can be constructed for each monthly data set using the technique 
of crossover constraints so that the variation of height at a common point in different solu- 
tions has an rms of around ±40cm. Some of this variation is probably due to oceanographic 
causes. 

There is no problem in generating orbits such that crossover discrepancies are not much 
greater than ±3m. These discrepancies are partly due to force model errors used in orbit in- 
tegration. They can be reduced to less than ±2m by introducing crossover constraints, pref- 
erably in terms of understood probable errors in the force field models, after elimination of 
about \0% of the altimetry data affected by excessive biases probably due to time tag errors. 
Given orbits which have a radial uncertainty of ±Ccm, further improvement can be achieved 
by averaging. The global data bank of sea surface heights (f or fr') can provide information 
for the recovery of a feature of the SST under the following conditions: 

(i) Sufficient GEOS-3 data is available for the average representation of any particu- 
lar wavelength so that the resolution is ±C/\/ncm, where n is the number of 
samples. 

(ii) The amplitude of the feature is greater than fi/y/n cm. 

(iii) The error in the geoidal model with wavelengths comparable to that of the feature, 
are significantly less than ±¥/v/ncm. 

It is not difficult to conjure up a scenario in which features of the SST can be deter- 
mined with a precision of ±10cm. For example, it has been observed that the second degree 
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of the SST is significantly larger than other terms in its representation by a surface harmonic 
series (Mather 1975, p. 67). On the basis of the previous study (Mather et al. 1977), it can 
be concluded that a basis exists for the recovery of this term from the GEOS-3 data bank, 
despite its flaws due to inadequate tracking distribution and the irregular manner in which 
the data were collected in both space and time. 

Other circumstances in which the magnitude of f 5 exceeds the noise level of the altim- 
eter are in the neighborhood of fast-flowing mid-latitude currents like the Gulf Stream. Whde 
there is little problem in obtaining f or f' on a relative basis in such circumstances with the 
required precision, the same cannot be said for f s due to the lack of knowledge of the geoid 
profile over such short wavelengths (less than 1000km), as discussed at length in (Mather 
et al. 1978b, Sec. 9). 

The present study therefore concentrates on the definition of parameters of the global 
quasi-stationary dynamic SST through those wavelengths for which the Earth's gravity field 
model is known with an acceptable level of precision. Section 2 deals with the data available 
for the present series of investigations. Section 3 deals with the technique for obtaining a 
solution from an analysis of long arcs of GEOS-3 altimetry. Section 4 deals with the Equinox 
Experiment designed to eliminate seasonal effects on the estimation of the quasi-stationary 
constituent of SST. The reference system plays a significant role in the evaluation of SST 
from a data set which is not truly global under circumstances where the solution can only be 
obtained for selected harmonics of the jp-avity field. This is discussed in greater depth in 
Section 5. Section 6 deals with the problem of modelling data which is restricted to ocean 
areas alone, with special reference to the use of the model in generating parameters related 
to ocean circulation. 

Section 7 analyzes the results obtained and Section 8 presents the conclusions drawn 
from the present study. 

2. THE 1977 GEOS-3 ALTIMETRY DATA BANK 

The orbits used in reducing GEOS-3 altimeter data were computed at Wallops Flight 
Center (WFC) using the GEM 7 gravity field model (Wagner et al. 1977) or at the Naval Sur- 
face Weapons Cen*e Dahlgren, Va., using tracking from 16 to 60 Doppler stations and a 
gravity field model specially tailored to the GEOS-3 spacecraft (Stanley 1978). Some orbits 
were also computed at Goddard Space f light C enter. The precision of the orbits provided 
by WFC appears to be better in 1976 than in 1975 (See Table 2). However, it cannot be in- 
ferred from these results that the precision of orbit determination is necessarily a function 
of time (ibid). 

In earlier studies (Mather et al. 1977, p. 34; Mather 1977. p. 24), it was shown that the 
application of corrections for a tilt c and bias b to each of a pair of “overlapping" passes over 
3000km lengths, made the median disagreement between passes drop to about t30cm. This 
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Table 2 

The Analysis of Crossovers (XO) for Twelve Ten Day Arcs used in the Equinox Experiment 

Wallops Orbits (See Section 4) 


Crossover Statistics 


Duration 


No. of 
Altimeter 
Passes 


Unadjusted 


No . of 
XO's 


Mean Root Mean 

L'iscrepancy Square Residual 


No. of 

XO's 


9. I .75 9.10.75 

9. 11.75 9.20.75 

9.21.75 9.30.75 

10. 1 .75 10.10.75 

10.11.75 10.20.75 

10.21.75 10.31.75 


+•3.49 


3. 1.76 


3.10.76 


Mean 

Discrepancy 

M 

+0.35 


+ 0.66 


3.11.76 3.20.76 52 53 +0.10 2.4 53 -°- 07 

3.21.76 3.31.76 46 28 +0.10 2.4 27 -0-35 

4. 1.76 4.10.76 34 30 +0.25 1.5 30 - O ^ 1 

4.11.76 4.20.76 29 10 -0.62 2.7 10 0.00 


12 4.21.76 4.30.76 j 40 63 -1.15 


Using the Orbit Error Model Defined in Equation 9 and a *5m cut oft w hich exduoed all passes aritb average XO discrepancies in excess of this limit 

















figure could be improved further if the slight non-coincidence of ground tracks were also 
allowed for (Mather et al. 1978b, Sec. 5). 

In the search for the dominant features of the dynamic SST, the primary objective is 
determination of the geometry of the sea surface in Earth space as defined in relation to a 
geocentric coordinate system. Such a model is directly compare ,c tothe geoid in any de- 
termination of f t . This would also apply to any satellite determined gravity field used in 
Ueu of the Utter in the search for dominant parameters defining the global distribution of 
{■,. This relation is conceptually d cribed by Figure 2. 

In view of the variable quality of the Wallops orbits, it was necessary to attempt an 
orbit improvement so that the radial orbital errors could be reduced 

(a) from values in excess of ±10m (e.g.. Mather vt al. 1977, p. 30) to values smaller 
than 1 2m; and 

(b) from values less than ±2 m to the relative precision of ±30cm implied in the anal- 
ysis of overlapping passes (Mather 1977, p. 25). 

The principal check on the quality of orbit integration is the crossover (XO) discrep- 
ancy d. While it is possible to achieve a median value of ±2m for d from orbit integration 
(Table 2), the improvement at (b) cannot be obtained from the GEOS-3 tracking data alone. 
The best tracking data available at the present time is laser data taken from the network of 
stations shown in Figure 3. The laser ranging precision varies from around ± 10cm for the 
NASA lasers to around ± \ m for other systems. The resulting orbit obtained from integration 
in the form of the instantaneous satellite position X, S (D at time t on a geocentric Earth space 
coordinate system X t , with the X, axis passing through the OO pole and the X t X, plane 
being that of zero longitude (X). The instantaneous position X iss (t) of the sea surface at the 
point with geodetic coordinates (0g.X) is obtained from the altimeter range h(t) using the 
relations 


X us (t) = X u (t)-C, h(t) + o{l0-*h} (4). 

where ^ are defined by the equations 

Cj = eosdgCosX Ci - eosdgSinX C* = sin (5). 

Timing is critical in correlating the orbital ephemeris with the time tag in the altimeter 
sensor data record as dh/dt can be as large as 20 m s' 1 . For example, a constant timing bias 
of flms causes sea surface height errors of 2d cm forcing soulh-to-north passes to have an 
error which is equal and opposite to that in a north-to-south pass at the same location. This 
has serious implications when attempting to enforce crossover constraints. 
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Equation 4 assumes that the altimeter is gravitationally stabilized. Orbit integration is 
usually performed over 5 day arcs using a gravity field model with the quality of Goddard 
Earth Model (GEM 9) (Lerch et al. 1977). The quality of orbit determination can be par> 
tially assessed by looking at the rms residuals for each span of tracking. While these are 
sometimes not as small as the quality of the tracking would warrant, this is largely attributed 
to along-track errors. This can be further verified by studying the quality of altimeter cross- 
overs for a 5 day span of data acquisition. A typical case is illustrated in Figure 4. 

The geocentric radial distance to the sea surface R s (t) is obtained from the coordinates 
X iss (t) using the relation 


R s (t) = 



( 6 ). 


R^ as evaluated at each altimeter data point during a 5 day arc, can be examined at crossover 
points for the propagation of the radial component of the orbital error as a function of time, 
using observation equations of the form 


R s (t,) - R s (t, ) + £c, 
1*1 


F *^2) - F »(t | > 


= V 


(7). 


where Cj are the coefficients required in the solution and Fj are functions of time. Observa- 
tion equations of this form could, in theory, include the ocean tide. In practice, however, 
the number of observation equations is rather small (less than 40 for a 5 day arc of GEOS-3 
data), precluding this possibility. It is important that crossover constraints be used only to 
eliminate orbital errors due to unmodelled force field effects. These are estimated as having 
predominant periods of one half revolution, one revolution. 14 revolutions, a resonance ef- 
fect with period of approximately 4.70 days and any linear drifts with time. 

Twelve ten day arcs of Wallops data selected for the Equinox Experiment (Section 4) 
were analyzed using crossover constraints per ten day arc using the following model for the 
orbital error 


) , CjF|(t) = C, t + C 2 cos 4 t r ™ + C, sin 4a ~ + C 4 cos 2*^- + C 5 sin 2tr— + 

j = j " t„ t„ t„ 


,, 2it t . „ . 2n t . „ 2tr t . „ 2n t , ox 

C* cos — 7 - + C n sm — — + C„ cos -t-tzt. — + C 0 cos — r. — ( 8 ). 


14 t_ 


14 t 0 8 


66.21 t 0 * 66.21 t 0 


where t is the time in days from the start of the 10 day arc and t 0 the orbital period in the 


10 







same units. The results obtained per 10 day arc are set out in Table 2. These results show 
that the root mean square (rms) residual of crossover discrepancies per 1 0 day arc averaged 
*6.8 m. This figure can be reduced to ±6.1 m by using equation 8 to represent the unmod- 
eled orbital errors (i.e., 18 percent of the power). A study of Table 2 shows that this average 
crossover discrepancy can be reduced significantly to ±1.8 by rejecting data associated with 
9 percent of the XO points on enforcing the criterion that any pass subject to an average XO 
discrepancy in excess of ±5m was subject to an unknown source of error, tentatively asso- 
ciated with time tag problems. 

Data sets referred to as “constrained” in this context are based on Wallops orbits as 
amended by XO constraints using equation 8 and the ±5m cut-off limit explained above. 

The significance of the coefficients Cj obtained by the use of equation 8 is not apparent at 
this stage. 

It therefore appears doubtful whether there is any means of obtaining profiles of sea 
surface heights from the orbital ephemeris of GEOS-3 using equation 4 with an absolute 
radial precision of better than ± 1.3 m in global terms at present. However, the internal pre- 
cision within the pass can be shown to be around ±30cm from overlapping pass comparisons, 
as discussed earlier. 

Data of this type can play a role in determining parameters of the global dynamic SST 
under the following conditions: 

(a) The error of ±1.3 metres mentioned above and which can be modelled by two 
parameters (a bias and a tilt) is randomly distributed as a function of position. 

(b) The parameter sought should not be much less than ±15cm. It should have a suf- 
ficiently long wavelength to enable its recovery from a global bank of data where 
satisfactory solution procedures can be devised for recovery of the dynamic SST 
signal under conditions where the signal to noise ratio approaches 0.1. 

If the sea surface topography were modelled using equation 43, the analysis of the data 
used to construct Figure 1 and summarized in Table 4 (Solution 3) shows that only the co- 
efficients ^ in and f sl20 satisfy the above criterion. However, if the orbital errors could in 
some way be brought to below ±60 cm there is some chance that the coefficients £ ri30 
and f d40 can also be recovered. 

Note that these probabilities are assessed only on the basis of an analysis of the GEOS-3 
altimeter data bank and do not take any other factors into consideration. This will be dis- 
cussed further in Sections 3 and 4. 
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Table 5 

I stimation of the Five Dominant Coefficients of the Quasi-Stationery Dynamic Sea Surface 
Topopraphy from (iEOS-3 Altimetry - Values (cm) 


Coefficient 


Oceanographic** 

(from 16,16) Solution 


Long Arc Solutions 


Constrained Wallops Orbits 


(_ iobal 


( 1 ) 

+114.5 
*6.9 
-21 .8 
-46.2 
*6.7 


.5 


Equinox 

Representation 

( 2 ) 

♦117.1 

♦9.6 

- 20.6 

-43.6 

+7.0 


Equinox Set 


1976 Data Set (till 8.76) 


AGS 638 
(3) 


EGM 19 
(4) 


Unadjusted 

(5) 

-27 

-153*** 

♦18*** 

-50* 

♦9 




Constrainedr 

( 6 ) 

-8 

-136*** 

♦ 6 *’* 

-35* 

-1 


Unadjusted 

(7) 

♦11 

-165*-* 

♦30*** 

-44* 

♦ 11 




Constrair^d T 

( 8 ) 

-4 

-141*** 

- 1 *— 

-35* 

-11 




rms voriotion 
per 1 °x 1° sq . 

± 5 

Average X0 
Discrepancy rms 
per 10 day ore 


No . of 1 ° blocks 

~ 




Percentage of Data 
Rejected (XO's) 


’These ..lues should be corrected by *6 1 tor me permanent t arth tide (Matner 19781 
**Bet««een 65 S and 65 N 

** 'includes non geocentrioTy o* the reference cooromate system 
t Using Equation 8 and '5m average XO limit 

























3. DETERMINATIONS BASED ON LONG PASS SOLUTIONS 


In attempting to reduce the average high level of orbital noise per pass from ± 1 .3 m to 
±30cm by introducing corrections for tilt c and bias b, the obvious goal is a means of deter- 
mining the corrections c and b without losing any information about the geometry of the 
sea surface. The introduction of the corrections c and b per pass effectively removes all in- 
formation with wavelengths greater than twice the length of the pass unless some special 
conditions were imposed. 

The longest passes of GEOS-3 altimetry are not greater than 9000km. Information on 
the shape of the sea surface which could be lost in tilting and bias correcting such a pass will 
be of wavelength 18,000km (i.e., harmonics less than degree 3 in a spherical harmonic rep- 
resentation of the SST). However, any contribution to the sea surface topography which is 
of degree 2 and symmetrical about the equator will not be lost on introducing corrections 
for tilt and bias for the following reasons. 

Let the ellipsoid of revolution which best fits the geoid have an equatorial radius a and 
flattening f. Let the difference in flattening between the ellipsoid of revolution which best 
fits the sea surface and the former be df. If the two ellipsoids have a common tangent at the 
equator (i.e., they have the same equatorial radius), the changes dR in radial distance R, 
given by 


R « a( 1 — f sin 2 # + ojf 2 }) (9), 

at latitude # between the two ellipsoids is 

dR « -a df sin 2 # + ojaf df} (10). 

dR is equivalent to f s if both the geoid and the sea surface had only second degree even 
zonal characteristics. The change SdR in the radial distance between two latitudes #j and 
<t> 2 is given by 


6dR = -adf(sin 2 # 2 - sin 2 #, ) (11), 

f„ 2 o being directly related to the second degree zonal harmonic in the pseudo-representation 
of the sea surface as a level surface. 

As * -46.2 cm (Table 4, Solution 3) the following conclusions can be drawn from 
a study of Table 3, noting that (e.g., Heiskanen and Moritz, 1967, p. 78) 

£ S 120 = df 02 ) ; 
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(i) No information on the meridian etlipticity of the sea surface or, for that matter, 
any other characteristic which is symmetrical about the equator (i.e., all even de- 
gree zonals in is lost by making corrections for bias and tilt to passes which 
are nearly symmetrical about the equator. 

(U) The lack of symmetry about the equator can be as large as 1000km for pa i of 
length greater than 6000km without causing more than a 50 percent distortion in 
the second degree harmonic. This is not considered acceptable for this type of 
work. However, unless some allowance of this type is made, it is not possible to 
sample the Indian Ocean region with passes of sufficiently long length. 


Table 3 

Differential Changes (Unsigned) in Ellipsoidal Height Due to Variation in 
Flattening over Finite Meridional Arc Lengths (s) (Equation 1 1)* 


Arc Lengths) 

Latitude 

6 dR 
(cm) 

Maximum 

6dR 

(cm) 

Minutes 

I'm 

Start 

End 

iO 

4425 

60 

20.2 

- .. 

96 

114 

40 

0.2 

63 

63 

20 

-19.8 

0 

18 

15 

6650 

60 

0.4 

114 

114 

40 

-19.7 

. - 

46 

63 

35 

-24.7 

24 

50 

20 

-39.7 

44 

62 

20 

8850 

60 

-19.5 

97 

114 

40 

-39.5 

1 

63 

39.8 

-39.8 

0 

62 

20 

-59.5 

95 

113 


"Computed ulinfl df - 2.381 * 10 7 (* •., t l)J0 " -46,2em) 
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Table 4 


Surface Spherical Harmonic Analysis of Oceanographically Determined 
Quasi-Stationary Dynamic Sea Surface Topography to (5,5' from Data Restricted 
to Oceans Between Parallels 65°S and 65°N (Units cm) 


Coefficients Normciiied 

Solution 

0) 

(2) 

(3) 

(4) 

Type 

Unconstrained 

Quadratures 

Constrained 

Quadratures 



Least Squares 

Zero on 

Least Squares 

Using Land 

n 

m 

Ocean Only 

land 

Ocean Only 

Bridge* 

— 

tt 

1 

2 

1 

2 

— 

1 

2 

1 

2 











0 

D 

110.2 


84.3 


114.5 


113.5 


n 

0 

5.3 


-8,5 


6.9 


3.4 


n 

i 

-27.9 

-2.7 

-34,0 

-4.1 

-21 .8 

+2.5 

-20.6 

+3.5 

2 

0 

-43.7 

wmm, 

-39.8 


• 46.2 




2 

1 

-13.1 

M 

-8.9 

-2,6 

-4.0 

+4.4 


+4.6 

2 

2 

40.8 

-4.6 

+5.7 

+1.0 

*0.7 

-0.2 

-0.1 

+0.1 

3 

0 

+ 18.5 


+6.5 


+6.7 


+ 1 .8 


3 

» 

-13.2 

-11.6 

+4.2 

-7.5 

-4.1 

-5.3 

-2.0 

-3.6 

3 

2 

+1 .8 


+8.5 

-13,6 

-0.? 

-2.4 

-0.0 

-1,9 

3 

3 

-0.7 


-3.2 

-10.4 

-3.0 

+ 1 .4 

-2.3 

+ 1.7 

n 

IS 

2.6 




-9.5 


-11 .9 


mm 

1 

-4.4 

+2.4 


1 3 .6 

+2.1 

+2.8 

+ 1 .9 

+2.7 

mm 


+2.1 

-6,8 


-0.9 


+ 1.2 

+ 2.7 

-0.6 

in 

WM 

+4.6 

-3.3 


+0.7 

+ 1.2 

-0.2 

+ 1 .5 

-0.4 

n 

Kfl 

-0.7 

-2,1 


-15.3 

-1 .8 

-0.9 

-1 .9 

-0.4 

5 


+ 14.7 


+ 2.7 


+ 1 .0 


-3.2 


5 

i 

-7.6 

+3.1 


+5.5 

-3.7 

-0.8 

-1 .0 

+0.4 

5 

2 

+1 .9 

-3.0 

+ 2.5 

+7.9 


+2.4 

-1 .2 

+2.6 

5 

3 

+4 .8 

-2.5 

-4.0 

+0.1 

+0.6 

+0.2 

+0.8 

+0.1 

5 

4 

+ 1.7 

+0.4 

-14.3 

+4.9 

-0.3 

+ 1.1 

-0.4 

+ 1 .8 

5 

5 

+0.6 

-0.1 

>0.8 

-6.3 

-0.3 

-0.8 

-0.4 


Anolysis 

To 

8,8 

16. 

16 




Rms fit 
(*em) 
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A basis therefore exists for recovering the dynamic SST from OIOS- 3 altimetry from 


long passes (lengths greater than s) under the following conditions: 
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(i) There should be sufficient passes of the maximum possible length, symmetrical 
about the equator, for the recovery of even degree constituents of the sea surface 
shape. 

(ii) Sufficient passes with a tolerance of about 500km from equatorial symmetry 
should be included to provide an even coverage of the Indian Ocean region. 

(iii) As will be seen from the discussion in Section 6, it is desirable to obtain a cover- 
age of the entire oceans to avoid aliasing effects due to incomplete sampling. 

Under such circumstances, there are several means of obtaining the necessary set of tilt 
and bias corrections. The simplest means, as argued above, is by fitting the set of long arcs 
to the best available geoid model, arc by arc, ensuring that the three conditions set out above 
are enforced. The following problems arise when using this method: 

(a) The resulting sea surface model will not have information on harmonies which are 
not symmetrical about the equator and with wavelengths greater than twice the 
length of the pass. 

(b) The harmonics of the SST deduced from a global set of such long passes of 

G EOS-3 altimetry are limited to those terms for which the error in the gravity 
field model is below at least 1 part in 10 8 . 

A study of this problem (Mather 1978b, Table 1) shows that, apart from considerations 
of the quality of the G EOS-3 altimetry data discussed above, the magnitude of the SST co- 
efficients and the estimated errors in the best available gravity field model - GEM 9 - indi- 
cate that only the coefficients f sl20 . f slll . f s!40 , f, |10 . f sl30 . f il60 and possibly, f sl21 and 
f s 22 i can be recovered at the present time. Of these, f s|11 . f sll() , f sI2 , and f s22 , cannot be 
recovered by the application of tilt and bias corrections. It is not clear whether a basis 
exists for the recovery of f sl30 from such data. 


Two types of solutions were obtained using the method of long passes which had 
orbital errors reduced from the ±1.5 m level to the +80 cm level by determining corrections 
b lor bias and e for tilt in relation to a gravity field model (GEM). The above three condi- 
tions for altimeter passes, although desirable to reduce correlation, are not necessary in this 
analysis since the bias and tilt parameters can be solved simultaneously with the SST coeffi- 
cients as in the solution AGS (>58 described later. These corrections were determined by 
setting up observation equations of the form 

v = h - t\ ;FM -|b + c(t-t 0 ) + 5h| (13), 

where h is the sea surface height obtained from GEOS- 3 altimetry, h GFM the height anom- 
aly at the sea surface computed by the application of Bruns’ formula to the GEM model. 



(t,t 0 ) being parameters quantifying the tilt in terms of height increments and 6h is the de- 
parture of the sea surface from the datum level surface (geoid), being modeled by the relation 

» ' in £ (r)" £ £ *■..« (14)> 

^o* ' O' <r»l 

where AC ftnm are coefficients of the surface spherical harmonic functions S Qnm defined by 
equation 27 and assumed normalized. The coefficients AC anm obtained from the resulting 
solution represent the geometrical distortion of the average sea surface from the level surface 
implicit in the GEM model. 

Equation 13, written in matrix form, is 

v » Ax + Bg- v 0 (15), 

where v 0 is the vector of differences (h - ), x being the vector of bias and tilt param- 

eters and g the array of coefficients AC anm . The reduced set of normal equations obtained 
from equation 15 by least squares, can be written in matrix form as 

[B T B-B T A(A T A)* l A T B]g * [b t - B T A(A T Af l A T ]v D (16). 

It can be seen from equation 13 that A T A for each pass of altimetry, is a 2 x 2 matrix 
with a simple inverse. Hence, for each pass(p) of altimetry, equation 16 may be written as 

V=R P (17) 

and summed over all passes to give 

Ng * R (18), 

If f is the matrix of weight coefficients which allows for variable observational accuracy 
and if solution constraints are introduced by the relation 


W g ~ 0 


(19). 


the final system of equations in matrix form are 

(f N + W)g - fR (19). 

The elements of the diagonal constraint matrix W (w fi corresponding to the element gj in the 
g array which is the coefficient AC anm in the representation of the SST (see Table 4)), are 
defined by the relation 
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( 20 ), 


Wj 


O * 
nm 


where 



o 2 being the degree variances (in cm 2 ) for the power spectrum of the dynamic SST. 


( 21 ), 


Two sets of solutions were attempted using this technique. 


Solution 1 (AGS 658) 

Data was selected for this solution to produce the best possible areal coverage. The 
distribution of data used is shown in Figure 10. 424 passes of G EOS-3 altimeter data, with 
ground tracks varying in length from 20° to 80°, comprised the data set. The input data 
were values of the sea surface height computed by Wallops Flight Center. Values of AC a(im 
were solved for to (5,5) using the following values for o n for purposes of constraint: 

o 2 = 1500cm 2 ; o 2 = 343cm 2 ; oj = 144cm 2 ; * 99cm 2 (22). 

For reasons put forward earlier, the coefficients AC^j and AC^ needed excessive con- 
straining, as did AC I33 , the respective values of w^ being increased by a factor of 3 over the 
values given in equation 22. The values used for f in equation 19 was 1 m. 

The GFM model used in AGS 658 (i.e.. for the generation of h CEM ) was GEM 7 
(Wagner et al. 1977) while values of the coefficients obtained in Table 5 are based on GEM 9. 
The degree variances in the uncertainties of the (5,5) model generated in solution AGS 658 
in comparison to those for GEM 9 are the following, in cm 2 : 


Degree 

AGS 658 

GEM 9 


66 

15 

3 

30 

125 

4 

'n 

66 

3 

12 

704 


On comparing the figures at (23) with those at (22), it can be concluded that starting 
with altimetry data with the level of noise around ±1 m, the technique has the potential to 
provide significant information on the quasi-stat ionary dynamic SST, subject to the limita- 
tions discussed above. However, the uncertainties in the GEM 9 coefficients, except for 
those of degree 2 and the zonal terms, restrict the information that tan be obtained at the 
present time. 
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Table 5 

Estimation of the Five Dominant Coefficients of the Quasi-Station&ry Dynamic Sea Surface 
Topography from GEOS-3 Altimetry - Values (cm) 


Coefficient 

Oceanographic** 
(from 16,16) Solution 

Long Arc Solutions 



Constrained Wallops Obits 

Equinox Set 

1976 Data Set (till 8.76) 

<_ ioba! 

Equinox 

Representation 

AGS 658 

EGM 19 

Unadjusted 

constrained* 

. J 

Unadjusted 

Constrained* 

n 

m 

0 

(1) 

(2) 

f3) 

(4) 

(5) | (6) 

(7) 

(8) 

3 

0 

1 

+ 1U.5 

+117.1 

- 

— 

-27 

-8 

+11 

-4 

i 

0 

1 

+6.9 

+9.6 

- 

- 

-153*** 

-136*** 

-165*** 

-141*** 

i 

1 

l 

-21.8 

-20.6 

- 

— 

+18*** 

+6*** 

+30*** 

*** 


0 

1 

-46.2 

-43.6 

-45* 

-47* | 

-50* 

-35* 

-44* 

-35* 

3 

0 

1 

+6.7 

+7.0 

+2 

+12 

+9 

-1 

+11 

-11 

- 

0 

s 

-9.5 

-9.9 

-7 

-26 

-32 

-38 

-29 

-28 

rnn vocation 
per !**1*jq . 

dfcS 

— 

±80 

— 

±172 

±149 

±184 

±161 

Average XO 
Discrepancy rms 
per 1C day ore 

_ 

— 

— 

— 

±680 

±180 

±332 

±177 

No. of 1° blocks 

- 

— 

— 

- 

13,499 

13,275 

12,349 

12,074 

No. of Data Points 

- 

— 

27, 67 4 

1757 

268,023 

232,687 

286,197 

234,929 

Percentage of Data 
Rejected (X0*s) 

_ 

_ 

— 

- 

- 

9 

— 

13 


■- — * — — - a ■■■Jtn — 

"These values should be corrected by *6.1 lor the permanent Earth tide (Mather 1978) 
“Between 65 S and 65 * 


** "includes non ^eocenir city of the reference coordinate system 
lUting Equation 8 and *5m average XO limit 




Solution 2 (EGM Series) 


In this type of solution, an attempt was made to recover primarily f j|20 using the tech- 
nique described by equations 13 to 21 . but with a smaller selection of passes restricted to 
those which satisfied the more rigid criteria set out following equations 9 to 1 2. One major 
problem was obtaining sufficient data to cover the Indian Ocean. In addition, the data sam- 
ple that satisfied the criteria of symmetry was about 40 times smaller than that used in 
AGS 658. This substantially increased the level of noise causing solution instability. 

Two types of solutions gave satisfactory values of f jiw . The first type (Solution EGM 
19) were passes longer than 4500km with balance across the equator to within 1500km. As 
only 21 passes of data in the 1977 GEOS-3 altimetry data bank satisfied this criterion, it was 
decided to assemble a second data set (EGM 20) which provided-, in addition a coverage of 
shorter passes (lengths greater than 1500km) in the Indian Ocean and Atlantic Region. The 
best results obtained from these two limited data sets were in the case where the higher ref- 
erence model used was GEM 9 to (30.30) where only zonal coefficients were solved for using 
a system of constraints defined by equations 19 to 22. The values obtained for f ji:o were 
-47 cm in the case of EGM 19 and -45 cm for EGM 20. Both these results are in 
agreement with the value expected from oceanographic considerations, despite the 
loss of resolution due to the limited data sample which is not well distributed globally (Fig- 
ure 8). Nevertheless the results substantiate the claim that information on the differential 
flattening between the sea surface and the geoid is not lost on application of corrections for 
tilt so long as the pass is symmetrical about the equatoi . 

An effort was made to improve the data coverage by reducing the criterion for equa- 
torial symmetry to 1000km. The resolution obtained deteriorated further, probably due to 
the aliasing effect of the additional number of passes asymmetrical about the equator. The 
result obtained from EGM 19 is shown in fable 5. though the values are not considered re- 
liable enough an estimate of the oceanographic value which is vulnerable to aliasing when 
evaluated from nor-global samples. 

This type of solution tor cannot be recommended due to the poor signal to noise 
and the paucity of data. However, they were made in an attempt to verify the thesis illus- 
trated in Table 3. that fbe application of tilt and bias corrections to passes symmetrical about 
the equator does not result in the loss of the signature of the second degree harmonic of the 
dynamic SST in the resulting sea s ivlace model. It can be concluded that this has been 
established. 

In conclusion, it can be stated that solution AGS o58 should provide information on 
the dominant harmonics i tXSr (possibly) and f s(4U of the dynamic SST. having a satis- 
factory areal distribution (Figure 10) and being confined to long arcs. The results obtained 
arc set out in Table 5. Solutions for harmonics of degree l ca* only be attempted from sea 
surface models related to geocentric orbits with appropriately controlled levels of noise. As 
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seen in Section 2, this is possible by applying some constraints when selecting a data base 
referred to the GEOS-3 orbits, as described in Section 4. 

4. THE EQUINOX EXPERIMENT 

One objective of satellite altimeter missions is the recovery of data for the synoptic 
monitoring of ocean circulation. As pointed out in (Mather 1978b, Sec. 4), over 80 percent 
of the spectrum of SST is contained in about 10 spherical harmonic coefficients through 
terms with wavelengths longer than 3000km. Attention is confined to these long wave fea- 
tures defining the global dynamic SST in this study. In the past, the lack of adequate d”ta 
has forced oceanographers to treat the global dynamic SST as stationary in time. This situ- 
ation has changed in the past two decades. 

The variations in the global SST with time can be classified as follows: 

(i) Synoptic variations with periods ranging from 10 1 to 10 2 days, which are prob- 
ably wind induced. 

(ii) Seasonal variations which are largely induced by the transfer of energy between 
the atmosphere and the surface layer of the oceans as a function of the seasons. 

(iii) Long period variations which are beyond the scope of the present study. 

The question of seasonal variation has been investigated by Gill and Niiler (1973). A 
comprehensive analysis of observational data has been undertaken by Wyrtki in the Pacific 
Ocean (Wyrtki 1975; Wyrtki 1977) and globally by Levitus and Dort (1977). 

The figures presented in these studies show that the changes in the SST are a function 
of season, sea level tending to be at its lowest in northern mid-latitudes during January - 
February and highest in July to August. The magnitude oi ihe variations is a function of 
position, varying from approximately 40cm in the Kuroshio to as little as 4cm at Guam 
(Wyrtki 1975, p. 457). The occurrence of such variations is of considerable significance in 
establishing sampling techniques for the determination of quasi-stationary SST. It is obvious 
that determinations of the global dynamic SST from data collected during a northern sum- 
mer will give different results to that collected during a northern winter. The terms affected 
will probably be the zonal harmonics, three of which contribute to the dominant 80 percent 
of the spectrum of sea surface topography. 

A separate determination of the dominant terms in the global SST ”'as therefore at- 
tempted from data which excluded samples recorded during the summer and winter months. 
This set is called the Equinox Data Set as only the short pulse mode altimetry data collected 
during the months of September - October 1975 and March - April 1976 were considered in 
its compilation. It was originally intended to use only the twelve 10 day spans of tta which 
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fell within the four month period mentioned. The selected U EOS-3 altimeter data distribu- 
tion for the September equinox in 1975 is shown in Figure 5 while that for t lie March 1 976 
equinox is shown in Figure 6. Fach distribution on its own is inadequate for determination 
of parameters of the global dynamic SS l . However the combination of the two (Figure 7), 
when supplemented by four additional five day spans on either side of the two periods, pro- 
vided a coverage of 39.8 percent of the 33,902 1° x 1° equi-angular blocks within the par- 
allels 65°S and b5°N which constitute the "global ocean areas" for the definition of the 
geoid used in this series of studies (Mather et al. 1978a, Sec. 0). 

Attempts were made to prepare a set of sea surface heights using the best available or- 
bits so that no information on low degree contributions to the shape of the sea surface were 
lost as in the case of the techniques described in the previous section. As discussed in Section 
2, the less than optimum distribution of tracking stations and limitations in the force field 
models can be expected to produce errors in the radial component of orbital position and 
hence, the sea surface heights. The use of crossover constraints provides a means of esti- 
mating some of the unmodeled errors in orbit integration, as seems to be the case from a 
study of the results in Table 2. 

The equinox data set was used both in the original and the constrained form, to prepare 
models of the sea surface as a first stage m the analysis lor those’ long wave components of 
the dynamic SST which have a favorable signal-to-noise in relation to the errors in the C'.FM 9 
gravity field model. The coefficients sought are set out in Table 5. The results obtained 
from this analysis are discusst'd in Section 7. I he data derived from the original Wallops data 
has a higher level of noise to that obtained after the application o' the error model at equa- 
tion 8 in the analysis of crossovers. In both cases the level of noise is significantly larger (a 
factor of at least 2 ) than the data sets used in Section 3. Operating under such circumstances 
is a necessary price which has to be paid in an attempt to retain all information on dynamic 
SST with wavelengths equivalent to spherical harmonics of degree less than 3 which are not 
symmetrical about the equator. 

Due to this adverse signal-to-noise, sea surface heights obtained directly from orbit in- 
tegration without the application of arbitrary tilt and bias corrections, require pre-processing 
into a suitable form prior to analysis. The nature of this pre-processing is described in the 
next section. 

5. fill SYSTFM OF RFFFRFNCT US! 1) IN Till EQUINOX I AIM RIM1 NT 


The principal role of a system of reference is the removal of systematic effects in the 
data which can be eliminated prior to analysis, thereby reducing the signal-to-noise ratio. It 
is also essential that the modelling procedure used does not inadvertently damp out the sig- 
nal sought. The model used in reducing the satellite altimetry data is that described as the 
higher reference model (Mather 1974. pp. 90 et seq ). I he use of this system in practice is 
described in (Mather 1 978b. Sec. 4). 
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In the fL;t stage, the height anomaly on the higher reference system is computed 
using the relation 


where 




( 24 ), 



(25), 


where (R o ,0,X) are the geocentric spherical coordinates of the sea surface, G is the gravita- 
tional constant and M the mass of the Larth. C^ nm are obtained from the satellite determined 
harmonic coefficients C anm to degree n' from the relation 


C' 

anm 


= C 


anm 


+ dC 


anm 


(26). 


As described in ( Mather et al. 1 978a, Sec. 6. 1 ), the corrections dC Qnm are primarily due 
to the differential effect of the atmosphere on downward continuation of the geopotential. 
These corrections turn out to be negligible for harmonics less than degree 6. The exceptions 
are the coefficients C' 120 , C' )40 and C' 160 . The use of the height anomaly ? in relation to a 
reference ellipsoid (equatorial radius a. flattening 0 in equation 24, requires that C' 1() = 0 
if the value of f is consistent with the coefficient C |1() . If this were the case, on forcing the 
above reference ellipsoid to be part of the higher reference model, the corrections dC, 40 and 
dC |60 take specific values which are functions of f (ibid. Appendix). Full details of the com- 
putation of f m are given in the reference quoted. The terms S anm in equation 25 are surface 
spherical harmonic functions, given by 


S lnm = P nm (sin0)c° s mX; S 2nm = P nm (sin 0) sin mX (27), 

P„ m (sin 0) being the normalized associated Legendre function of degree n and order m. 

The basic relation satisfied by the height anomaly f' at the surface of the Farth is 
(Mather 1975, p. 72) 


T = ( w 0 -u 0 ) + 7 r-7i: s 


(28), 


7 being normal gravity. 

The difference (W t) - U c ) is established from the GLOS-3 altimetry as described in 
(Mather et al. 1978a, Sec. 6) and is a known quantity, T, for all practical purposes is the dis- 
turbing potential which will be zero for all harmonics less than n' if the coefficients C onn) 
are also free from error, the atmospheric potential being neglected in equation 28 as its dif- 
ferential effect has been shown to have no significant influence on relations of low degree at 
the surface of the Larth not involving upward continuation (ibid, p. 17). 
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A solution for low degree harmonic coefficients (n < n') in f s can only proceed on the 
basis that the only other term in equation 28 containing such information, is the global data 
bank of f' which, in the present investigation, is only available at sea between parallels 6S°S 
and 65° N. 

The input data in the Equinox Data Set is in the form of 13,499 1° x 1° area mean 
values of f' (Table 5), with a rms variability of ±2.7 m. The mean variability within a 1° x 1° 
square from raw Wallops orbits is ±1.7m. As the harmonics being sought in the present study 
are of degree less than 7, it is preferable to analyze the data in the form of 5° x 5° area means 
as the signal-to-noise is improved. 

Simplified observation equations of the form 

v = * - L £ H f sanm S «nm < J 9 ), 

n-1 m=0 o=l 


where f sanm are the surface spherical harmonic coefficients in the global representation of 
f s , k being defined by either 


k =-(w 0 -u 0 )±f' 


(30) 


or 


k = - 
7 


w ° ' K It, (0 tk !l <a,,m S “"'" ' 2 R » cos: * 


( 21 ). 


U Q being the potential of the rotating equipotential ellipsoid defined by the parameters GM, 
a, f and to, where u> is the angular velocity of rotation of the Earth by the relation (e.g.. see 
Mather 1971, p. 83) 


where 


U. 


GM 


a 


a 

sin a 



a 2 ur 


(32). 


a = cos* 1 ( 1 - f) 


(33). 


The satellite altimetry data is used to define the geocentric spherical coordinates 
(R o ,0,X) of the sea surface when using equation 31 to define k. 

The reference system used for the Equinox experiment was defined by the following 
set of constants: 
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( 34 ). 


c » 2.997 924 58 x 10 10 cm s’ 1 
GM = 3.986 004 7 x 10 :o cm 3 s’ 2 
C pu = 1.082 627 5 x 10° 
u> * = 7.292 115 15 x 10* 5 rads*' 

The dependent constants have the following values: 

r 1 * 298.257 3 (35). 

The potential of the geoid as obtained from the analysis of the equinox data (Mather et al. 
1978a. Sec. 6) using the constrained set from the Wallops orbits (Table 5), 

W 0 = 6,263,682.76 ±0. 1 8 kGalm (36). 


giving 


(W 0 “U 0 ) = +0.08 kGal m (37), 

for an ellipsoid of equatorial radius 6,378.140.00m. 

The geoid so defined is for epoch 1976.0. It should be noted that W 0 is based on a 
39.8 percent representation of the ocean areas between 65°S and 65°N. with data distribu- 
tion as shown in Figure 7. The calibration adopted for the GEOS-3 altimeter is that deter- 
mined pseudo-geometrically by Martin and Butler ( 1977). This procedure, while affecting 
the value of W o does not influence the dynamic SST parameters determined from the 
Gl OS-3 data bank. 

6. MODELLING fill DYNAMIC SIA SURF ACL TOPOGRAPHY 

Special problems are involved in modelling data which does not continuously cover the 
surface of the Earth. Oceanic phenomena cover only about 70 percent of the Earth’s sur- 
face. The data for GEOS-3 only provides data for those parts of the oceans which lie be- 
tween 65 *S and 65°N. On using the 1" \ I global elevation data bank as a mask for the 
oceans, this defines the ocean area in terms of 33.902 such squares between the parallels 
defined above. 

Extractable information on the quasi-stationary dynamic sea surface topography is 
contained in the GFOS-3 altimetry data bank subject to the following qualifications: 

(i) The features sought should preferably have a magnitude in excess of at least three 
times the noise level of equivalent wavelengths in the gravity field. 

(ii) Sufficient altimetry data is available to filter out the effect of the tides. 
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(iii) Adequate data is available to average out any seasonal variations in the SST. If 
this were not the case, the resolution obtained will only be equivalent to the mag- 
nitude of such variations (currently estimated at ± 20 cm). 

The obvious model for representing the long wave features of the SST is a surface 
spherical harmonic series 

2-! ^sanm ^onm + °{^s} (38), 

n*l m*0 o«l 

where f >anm are coefficients of degree n and order m, while S ftnm are defined by equation 
27. The values obtained for the coefficients f sanm are dependent on the method of solution. 
The harmonic coefficients can be obtained by two different methods if the data is uniformly 
distributed over a sphere. 

(i) By least squares, using observation equations of the form 

AX-K = V (39), 

from equation 29. Values of f sanm are obtained by minimizing 

q, = V T W V (40), 

where W is the matrix of weight coefficients which in this set of computations is 
cos 0, where 0 is the latitude of the equi-angular square in which the data was 
sampled. 

The solution is then 

X * (A T WA)-'A T WK (41), 

The array X being composed of the coefficients f sanin . 

(ii) Alternately, given a distribution of f s , the coefficients f sann , can be obtained 
using the relation 


^sanm 4 ^ J J ^anm (42), 

assuming the coefficients to be normalized, do being the element of surface area 
on unit sphere. 
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The error introduced by the spherical assumption is less than ±0.5 cm as the SST does 
not have a magnitude in excess of ±2ni. 

The question of modelling the long wave features of the SST has been discussed before 
(Mather 1975; Mather 1978a. Sec. 6). Two distinct possibilities are open: 

(a) Solve for f SQnm using equations 28 to 31 , sampling ? s in ocean areas only. 

(b) Solve for f SQnm using equation 29 but replacing by 

«■ 2 

f s - k o (0,X) X- ^sanm ^anm (43), 

n=l ni=0 a*l 

where 

k () (0.X) is the ocean function defined by 
f 0 if t^.X ) on land 

k o (0.X) = s (44). 

I 1 if (0.X ) is oceanic 

Table 4 sets out the coefficients obtained by analyzing the quasi-stationary dynamic 
SST represented in Figure 1 by four different methods. Solution 1 gives the results obtained 
by using equation 29 in ocean areas only. Solution 2 gives values of ? SQnm obtained by using 
equations 42. 43 and 44. The results for solution 2 change negligibly between equivalent 
least squares solutions and values obtained by quadratures. The results in solution 1 for a 
determination to (tS.8) vary considerably when the analysis is carried out by least squares to 
differing degrees, due to a high level of correlation between the solved coefficients as a re- 
sult of the incomplete representation of data on the sphere. 

A third type of solution was obtained by constraining the coefficients f SQnm bj im- 
posing the condition 


4* = V T WV + X T W C X = Minimum 


(45) 


instead of equation 40 when obtaining a solution. The constraints were imposed through the 
array W , using coefficients of the form 


W 


1111 



10 3 


(4b). 


for answers in in. where n is the degree of the harmonic and R the mean radius of the Earth. 
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Solutions obtained in this manner changed only slightly as a function of the degree to 
which the data was analyzed. A solution of this type is shown as solution d of Table 4. 

The constrained solution in ocean areas so obtained was checked by interpolation ol 
values in laud areas along a parallel using a cubic polynomial in longitude difference d\ from 
the eastern sea/land boundary, of the form 


•i 

a, <dX)‘ (47). 

1-0 

The coefficients a, are estimated by least squares fitting to between 4 and 8 data points 
in oceans along the same parallel on either side of the land mass which needs to be bridged. 

The results so obtained from both least squares and quadratures are listed as solution 4 
in Table 4. 

The results in Table 4 emphasize the care needed in selecting a system for modelling the 
SST. No two models in Table 4 are exactly equivalent , though the models in the last two sets 
of columns differ only in that zero has been assumed for values of f s outside the parallels 
b5 J N and b5°S in one evaluation by quadratuies. 

The ^.efficients in unconstrained solutions are heavily correlated and change dramat- 
ically with tl.e maximum degree to which analysis is carried out. The quadratures solutions 
using zero on land are not directly comparable with harmonic solutions obtained front sets 
of observation equations formed in ocean areas only, furthermore, the use ot such models 
in current modelling li quation 1 ) in the form of a horizontal derivatives, become unstable 
close to coastlines, as do the representations of due to forcing the harmonic model to 
take zero values on land. 

The variation in the coefficients for solutions of types .) and 4 indicate the strong in- 
fluence of areal representation in solutions, l i I OS-.) altimeter data docs not provide values 
of sea surface heights in all oceanic areas between b5°S and bV‘N. as can be seen from the 
distribution of data for the equinox experiment (figure 7). 

It is therefore considered essential that this factor be taken into account when assessing 
the relative quality of solutions obtained from different incomplete data banks. The domi- 
nant surface spherical harmonic coefficients for the quasi-stationary SST used in the present 
study tti represent the fully sampled region between the two sixty-fifth parallels in ocean 
areas is given in Column I of Table 5 Column - ol this table gives the values obtained had 
the analysis been done using data sampled only in the areas where altimetry data were avail- 
able for the I quinox experiment. Moth sets of values are obtained from constrained least 
squares solutions to ( In. lb). 
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7. HU ANALYSIS Ol III! KISl'l IS 


Table 5 presents in summary, the results obtained to date in determinations of the 
dominant parameters of the dy namic sea surfaee topography (SS L). The results show that 
correctly filtered Ci EOS-3 altimetry data can provide heights of the sea surfaee in relation to 
a geocentric coordinate sy stem with a sufficiently large signal to noise to permit the recovery 
of the parameters referred to above The level of noise in the constrained data sets obtained 
from the Wallops orbits using equation 8 on 10 day ares, is estimated from crossover (XO) 
discrepancies at ±l.3m. A second indicator is the variability of values in a 1° x I" square. 

The mis geoid variation across such a square varies between ±0.3m to iO.Hm. The value of 
±1.5 m obtained for data in the equinox data set seems to confirm the estimate of ±l.3m as 
the level of noise in the data comprising the I quinox Data Set 

It has been a rule of thumb to assume that it is possible to derive long wave signals 
whose amplitudes are ten times smaller than the noise level provided the background noise 
does not contain errors of equivalent wavelength and amplitude. The analysis of XO’s to 
model short period orbital error .in described in Section 2 , reduced the XO discrepancies on 
average by ±3m. The resulting constrained solution produces a value of f g|J0 which is sig- 
nificantly smaller than the oceanographic value. This casts a shallow over any ad hoc model- 
ling of terms in the enforcement of XO constraints. 

In the present set of solutions, a limit of * I m were set on the magnitude of the coeffi- 
cients C'j m equation S, using equation 45 I Ins may be too high Therefore, the constrained 
set of solutions m the present study may not necessarily be the best solutions. 

The analysis of Xl)*s for the entire l l >77 i ;| OS-3 altimeter data set in 10 day spans in- 
dicated that the orbits in l l >7o were significantlv better than those in 1975. It was therefore 
decided to compile a l°7o Data Set based on data ..equired between January and August 
l‘>7{» i Figure 9). The same two tv pes of solutions as obtained for the I quinox Data Set 
were repeated and the results listed m the last two columns of I able 5. \gam. the constrained 
solution gives a relatively low value of , 0 , confirming the need for a review of the proce- 
dure for modelling orbital error in enforcing XO constraints. 

I he Second Degree Harmonic 

In view of the reservations about the constrained solutions, the authors are inclined 
towards a value of -47cm for f However, this value is only a preliminary estimate. It 
should also be noted that the value is obtained by referring the true sea surface to a static 
I arth model which is independent ot the I artli tide The latter has a permanent component 
which can be estimated as contributing to the shape of the sea surface A value of v com- 
parable w ith the oceanographic values is obtained In correcting the above value by ♦ocm 
t Mather 1978c). 


33 




Figure 9. Distribution of GEOS-3 Altimetry in the 1976 Data Set 





e 10. Distribution of GEOS-3 Altimetry Used in Solution AGS 658 





The First Degree Harmonic 


It had been hoped to recover the second largest contribution to the dynamic SST 
(f NllI - Table 5) which has a value -21cm. Being a first degree harmonic term, it can only 
be obtained from the analysis of the orbit related sea surface height data bank (i.c., either 
the Equinox Set or the 1976 Data Set). Solutions for this harmonic can only be obtained if 
the GEOS-3 orbits refer to the geocenter without error. This is not the case. As the first 
degree harmonic of f s has a degree variance of 5 1 5 cm 2 , it follows that the center of the 
ellipsoid of best fit to the sea surface does not depart from the geocenter by more than 
25 cm. The degree variance of the value obtained from the Wallops orbits in all cases exceeds 
1 8500cm 2 . The only inference that can be drawn is that the origin of the system of refer- 
ence used in integrating the orbits is displaced from the geocenter by not more than ±1.5m. 

The largest contribution, by far comes from the first degree zonal harmonic (Table 5), 
which indicates a southward displacement of 1 .49 ±0.1 5 m. As the sources from which or- 
bital infom.ation were obtained varied, no specific conclusions are drawn from the^e num- 
bers at the present time. 

The Zonal Harmonics of Degree 3 and 4 

The magnitudes of f sno and $ sl4n are smaller than one-tenth the estimated systems 
noise in the case of the orbit related SST. The signal to noise approaches the 0.1 level in the 
case of the long arc solutions AGS 658 and EGM 19. As pointed out earlier, it is debatable 
whether the long arc solutions can provide estimates of f sl30 . Disregarding values from the 
suspect constrained solutions, the average result from Table 5 is 

f*i 3 o = ±5 cm (48). 

The standard deviation is significantly smaller than expected from the input data noise 
levels and is in very good agreement with the oceanographically determined value. If the 
standard deviations are to be taken at face value, the seasonal effect on f s , 30 should be absent 
m the result from the equinox experiment. The difference of f2cm is consistent with the 
1976 data set representing data during a southern summer. The irregular distribution of data 
between the hemispheres (Figure 9) cautions against any drawing of conclusions from this 
result. 

The value of ? sl4(l is not so well defined. The average value obtained without enforcing 

crossover constraints is 

f s)40 = -23 ± 1 8 cm (49). 

The standard deviation is much more in keeping with that estimated from the noise of 
the input data. The values obtained from orbital and long arc solutions differ somewhat. 

The reason for this difference is not clear at this stage. 
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8. CONCLUSIONS 


The following conclusions can be reached: 

(i) A basis exists for determining three out of the five dominant parameters of the 
quasi-stationary dynamic SST from GEOS-3 altimetry despite the fact that the 
data acquisition patterns did not have global oceanographic objectives in view. 

(ii) The present series of solutions are confined to the basic orbits generated for 
Wallops Flight Center, as modified in ways prescribed above. The results obtained 
are preliminary and subject to revision as improved orbits become available. On 
the basis of the discussion given above, the best values for the dominant harmon- 
ics other than those of degree one, in the SST are the following 

f sPn (Second Degree Zonal) = -43 ±6cm (-46) 

? S | 3 o (Third Degree Zonal) = +7±5cm ( +7) (50), 

f sl40 (Fourth Degree Zonal) = -23 ± 18cm (-10) 

the values in brackets being the oceanographic estimates. The value for f sP() has 
taken into account the contribution of the permanent Earth tide and the -2cm 
effect introduced by the irregular data sampling (Table 5. Columns 1 and 2). 

(iii) It is hoped to improve on these results by a revision of the constraining procedure 
used in preparing Table 5 (i.e., equation 8) in conjunction with improved laser 
supported orbits. 

(iv) An improved gravity field model for the low degree tcsseral terms aiming for a 
resolution of 3 parts in lO 4 is an important pre-requisite to further progress in 
this area. 

(v) It cannot be too strongly emphasized that an evenly spaced altimeter data cover- 
age of the world’s oceans is indispensable for recovering the long wave features of 
the quasi-stationary sea surface topography with confidence. On the basis of the 
present study, this could be achieved if the radial component of orbital precision 
can be improved by one order of magnitude. 
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